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ABSTRACT 

As Lyman a photons are scattered by neutral hydrogen, a change with redshift in the Lya equivalent 
width distribution of distant galaxies offers a promising probe of the degree of ionization in the 
intergalactic medium and hence when cosmic reionization ended. This simple test is complicated 
by the fact that Lya emission can also be affected by the evolving astrophysical details of the host 
galaxies. In the first paper in this series, we demonstrated both a luminosity and redshift dependent 
trend in the fraction of Lya emitters seen within color-selected 'Lyman-break' galaxies (LBGs) over 
the range 3 < z < 6; lower luminosity galaxies and those at higher redshift show an increased likelihood 
of strong emission. Here we present the results from much deeper 12.5 hour exposures with the Keck 
DEIMOS spectrograph focused primarily on LBGs at z ~ 6 which enable us to confirm the redshift 
dependence of line emission more robustly and to higher redshift than was hitherto possible. We find 
54 ± 11% of faint z ~ 6 Lyman break galaxies show strong (WL ya ,o > 25 A) emission, an increase of 
1.6x from a similar sample observed at z ~ 4. With a total sample of 74 z ~ 6 LBGs, we determine 
the luminosity-dependent Lya equivalent width distribution. Assuming continuity in these trends to 
the new population of z ~ 7 sources located with the Hubble WFC3/IR camera, we predict that unless 
the neutral fraction rises in the intervening 200 Myr, the success rate for spectroscopic confirmation 
using Lya emission should be high. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: starburst - galaxies: high 
redshift 



1. INTRODUCTION 

The reionization of neutral hydrogen in the intergalac- 
tic medium (IGM) was a landmark event in cosmic his- 
tory, rendering the Universe transparent to UV photons 
and dramatically reducing the star formation efficiency in 
dwarf galaxies. In spite of its importance, there are few 
robust constraints on when reionization occurred. Po- 
lariza tion measures of th e microwave background radia- 
tion (|Larson et alJ feOlO) demonstrate scattering by free 
electrons in the redshift range 7 < z < 20 but do not 
describe the evolving neutral fraction, Xhi- Absorption 
line spectra of high-z quasars are largely sensi tive to the 
very late stages (xhi — 10 -3 ) of reionization (|Fan et al.1 
|2006T >. and progress has been slow due to the paucity of 
sources so far detected beyond z ~ 6.5. 

One of the most promising probes of reionization with 
current facilities is through the study of Lya emission 
from star forming galaxies. Since Lya photons are reso- 
nantly scattered by neutral hydrogen, the abundance of 
Lya emitters should decrease as observations probe into 
the era where there are pockets of neutral gas. Stud- 
ies of the redshift-dependent luminosity function (LF) 
of Lya emitters (LAEs) selected via narrowband filters 
have revealed a possible decline in abundance between 
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z = 5.7 and z = 7.0 dKashikawa et al l [20061: live et al.l 
120061: lOta et~alll2008l: lOuchi et aJj|2oToji . offering tanta- 
lizing evidence that this short time interval (~200 Myr) 
may correspond to one during which there is some evo- 
lution in the neutral fraction. But since a number of 
astrophysical factors can also affect the presence of Lya 
emission, it may be dangerous to dir ectly link evolutio n 
in the Lya LF to reionization (e.g., iDaval et al.l [2010). 
These factors include time-dependent changes in the host 
galaxy number density, dust obscuration and interstellar 
gas content and kinematic properties. By enlarging the 
LAE samples, it may be possible to bypass some of these 
complications by testing for the expected change in their 
spatial c lustering and line profiles as the neutral era is 
entered (|Ouchi et al.ll2010D . 

A c omplementary approach introduced in IStark et al.l 
(2010) (hereafter Paper I) is to spectroscopically mea- 
sure the fraction of strong Lya emitters within the color- 
selected Lyman Break Galaxy (LBG) population. By 
tracing the redshift-dependent fraction, the host galaxy 
number density is not a factor. Evolution in dust ob- 
scuration can be independently tracked using the con- 
tinuum c olors and ISM kinematics through deep spec- 
troscopy (ISteidel et al.ll2010HBouwens et al.ll2009ll2010al: 
iVanzella et al.ll2009() . Although demanding observation- 
ally, high throughput spectrographs such as FORS2 on 
the ESO Very Large Telescope and DEIMOS on the Keck 
II telescope have enab led progress in recent years (Paper 
I, IVanzella et al.|[2~009f) . With the additional information 
on the host galaxies possible for the LBG population, we 
can hope to more reliably link any redshift-dependence 
in the Lya fraction to ionization changes in the IGM. 
Most importantly of all, the proposed approach can now 
be readily extended to z ~ 7-8 and beyond given the 
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availability of LBG samples at these early epochs follow- 
ing the advent of the WFC3 /IR camera onboard H ubble 
Space Telescope (HST, e.g. iBouwens et~alll2010al ). 

In Paper I we introduced a large Keck spectroscopic 
survey of z > 3 LBGs and demonstrated the practi- 
cal details of the above method through analyses of the 
Lya fraction (XLya) m B-band (z ~ 4) and V-band 
(z ~ 5) dropouts to which we added a sample of i'- 
band fz~6) dropou ts drawn from other programs (e.g., 
IVanzella et ai1l2009l Bunker et al. 2010, in prep). Cor- 
recting for minor magnitude and redshift-dependent bi- 
ases in completeness and contamination, we determined 
the luminosity and redshift dependence of Xl vol over 
3 < z < 6. Since the IGM is known to be highly 
ionized over this interval, this dataset enabled us to 
explore the importance of factors other than the IGM 
neutral fraction. We found that galaxies with lower 
rest-frame UV continuum luminosities exhibit Lya emis- 
sion more frequently than luminous systems. Correla- 
tions between line strength and UV continuum slope 
suggest reduced dust obscuration is the primary cause. 
The data also suggest an increase in XL ya with redshift 
(dxLya/dz ~ 0.05 ± 0.03), as originally claimed based on 
the relative evolution of the UV luminosit y function of 
LAE s and LBGs over this redshift range ()Ouchi et al.l 
2008). However, since the size of our archival z ~ 6 
samples were considerably smaller than that of our lower 
redshift database, the constraints on the redshift evolu- 
tion were primarily derived from data spanning only the 
300 Myr between z ~ 4 and z ~ 5. 

Ideally one would construct the full equivalent width 
(EW) distribution function of Lya emission for a range 
of UV luminosities at the highest redshift where the IGM 
is known to be highly ionized, i.e. z ~ 6. This could then 
form the basis for comparisons with spectroscopic data 
at z > 7 where reionization may be incomplete. With 
this motivation, we have thus extended the sample in- 
troduced in Paper I, through ultra-deep spectroscopy of 
a sample of i'-band dropouts in the GOODS North held. 
The increased sample size and deep spectroscopic expo- 
sures provide statistically-significant constraints on the 
EW distribution of feeble sources, ensuring an adequate 
basis for comparisons with higher redshift spectroscopic 
samples. 

Throughout the paper, we adopt a A-dominated, 
flat universe with f^A = 0.7, Hm = 0.3 and H = 
7OI170 kms _1 Mpc _1 . All magnitudes in this paper are 
quoted in the AB system (|Oke fc Gunnlll983l ). 

2. OBSERVATIONS 

Our dataset is primarily comprised of spectra ob- 
tained using the DEep Imaging Multi-Object Spectro- 
graph (DEIM OS) at the Nasmy th focus of the 10 m Keck 
II telescope (jFaber et al.l 120031) . We direct the reader 
to Paper I for a full description of our survey strat- 
egy. In Paper I we presented analysis of 513 DEIMOS 
spectra, including 268 unique B-drops and 95 unique 
V-drops. To this we added publicly available spectra 
from the VLT/FORS 2 survey of z ~ 4, 5, and 6 LBGs 
(|Vanzella et al.ll2009|) and 2 unique z ~ 6 LBGs from the 
Keck survey of Bunker et al (2010, in prep). The to- 
tal sample drawn from Paper I is thus 351 B-drops, 151 
V-drops, and 44 i'-drops. 

The major step forward here is the inclusion of new 
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Fig. 1. — Montage of 2-D Lya detections in the April 2010 
DEIMOS run targeting i'-band dropouts in the GOODS North 
field. The color scale is inverted with black corresponding to pos- 
itive flux. Each cutout spans 7.1 arcsec X 28 A. The i'-band 
dropouts that we observed are faint, with optical magnitudes span- 
ning 26.0 < 2850 < 27.5. The Lya detections are coincident with 
the spatial position of a UV continuum dropout satisfying the i'- 
band dropout color criteria (i'-z > 1.3 and containing no detec- 
tions in deep B435 or V606 imaging). Only emission lines with 
WLyafi > 25 A are included in the analysis in §3 and §4 (which 
excludes the second from right source in the bottom row) , ensuring 
that the EW distribution is derived from robust detections. 



z ~ 6 spectra following ultra-deep Keck exposures of 
faint i'-band dropouts in GOODS-North. The archival 
data in Paper I was mostly based on the equivalent of 
3-4 hour exposures with a 10 meter aperture. The new 
sample consists of 23 z'-band dropouts with 12.5 hour 
exposures (and an additional 7 with 3.67 hours of inte- 
gration) enabling constraints to be placed further down 
the EW distribution at z ~ 6 (allowing a uniform sam- 
pling over the full redshift range) and increasing the total 
z ~ 6 LBG sample by nearly ~ 70% to 74 across both 
GOODS fields. 

The new data were taken during April 2010. Over 
11-12 April, we obtained 12.5 hours of on-source inte- 
gration in good seeing (<0'.'8) for one mask containing 
23 i'-band dropouts. Over 13-14 April, we obtained 3.67 
hours of integration on a separate mask containing 7 i'- 
band dropouts. For both masks, we used the 830 line 
nun" grating, typically providing spectral coverage be- 
tween 7000 A and 10400 A. Slit lengths were gener- 
ally ~ 7", and slit widths were 1". Skylines are mea- 
sured to have a Gaussian a of 1.1 A. Reduction was per- 
formed using t he spec2d IDL pip eline developed for the 
DEEP2 survey (|Davis et al.ll2003l ) 7 . Wavelength calibra- 
tion was performed using Ne+Xe+Cd+Hg+Zn reference 
arc lamps. As in Paper I, we flux calibrate our data using 
the spectra of alignment stars included on the slitmask 
(observed in 2" by 2" boxes). We compared this cal- 
ibration to that obtained using spectroscopic standard 
stars and found it to be consistent to within ±20% (with 
no significant systematic offset) for the alignment stars. 
Using the flux calibration, we computed our survey sen- 

7 The spec2d pipeline can be downloaded at 
http : / /j uno.as . arizona/cooper /deep / spec2d / 
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Fig. 2. — Left: The differential rest-frame Lya EW distribution, p(WLya) (computed in bins spanning AWLya,o = 30 A) for star-forming 
galaxies at z ~ 6 in two luminosity ranges (—21.75 < Mjjv < —20.25 on bottom and —20.25 < Mjjv < —18.75 on top). Overplotted in 
red is the Lya EW distribution for LBGs at 4 < z < 5 derived from the sample in Paper I (dotted lines provide la uncertainties). Right: 
Evolution in the overall fraction of Lya emitters (Xi 9a ) in the LBG population over 4 < z < 6. Luminous LBGs are considered in the 
bottom panel, and less luminous systems in the top panel. In each panel, we derive the Lya fraction of LBGs with Lya EWs larger than 
25 A (circles) and 55 A (squares). Assuming a linear relationship between Xi Btl and z, we extrapolate these trends to z ~ 7 (triangles 
with dashed-line error bars). 



sitivity as a function of wavelength. The 5a limiting line 
flux is 3.1±0.5xl0 -18 erg cm -2 s _1 (assuming a range 
of Lya line widths typical of our LBG samples) , imply- 
ing that we should be able to detect Lya with rest-frame 
equivalent widths of greater than 20 ± 3 A for i'-drops 
with zg5o — 27. 

3. ANALYSIS 

We searched for Lya emission at the spatial position 
of the targeted LBGs in the Keck spectra. Line fluxes 
and EWs were calculated following the procedures dis- 
cussed in Paper I. We account for the effects of line con- 
tamination and Lya forest absorp tion (estimated using 
relations presented in|Meiksin 200§) on the observed z§5o- 
band fluxes. Of the 23 i'-band dropouts for which we ob- 
tained ultra-deep spectra, 11 show Lya emission, while 
2 of the 7 i'-drops for which we obtained 3.67 hour in- 
tegrations show Lya (Figure [1} . These results imply a 
large fraction of i'-drops have prominent Lya emission. 
The rest-frame EWs for the i'-drops range between 9.4 
A and 350 A. The vast majority of the emission lines 
are detected with high significance. Even so, we take a 
conservative cut, limiting our analysis to those sources 
with rest-frame EWs greater than 25 A and S/N > 7. 
This excises the Lya detection in the middle right bottom 



panel of Figure [T] As a result, even among the faintest 
sources, the emission lines used in our analysis are very 
confidently detected (< S/N >=18), removing concern 
regarding spurious features. 

As in Paper I, we determine the completeness of our 
Lya detections as a function of wavelength by adding and 
recovering fake emission lines at random positions across 
the 2-D spectra. We compute the Lya recovery rate as 
a function of absolute magnitude and wavelength for all 
masks observed (including those in Paper I) and make 
appropriate corrections. This test demonstrates that in 
our deep 12.5 hour mask we are >90% complete to lines 
with Wiy Qi o > 50 A even for the faintest i'-drops on 
our mask (z$ 50 ~ 27). For lines with WL ya ,o — 20 A, the 
completeness implied by our simulations is ~ 75 — 80% 
for sources in the faintest magnitude bin covered by our 
z ~ 6 spectra. The completeness is of course lower on 
the DEIMOS mask observed for only 3.67 hours, reach- 
ing below ~ 50% for faint sources with WL ycti o — 20 A, 
and we therefore do not include these sources when com- 
puting the fraction of LBGs with low EW Lya emission. 

An additional concern is that the color-cut and z-band 
selection of z'-band dropouts are affected by Lya emis- 
sion and Lya forest absorption. We investigate the ex- 
tent to which these effects transform the observed EW 
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distribution using Monte Carlo simulations. We create 
a large sample (>10 6 ) of artificial galaxies with intrin- 
sic absolute magnitudes (normalized at 1500 A) span- 
ning —21.5 < Mjjv < —18.5 and redshifts spanning 
5.6 < z < 6.5. The intrinsic luminosity distribution 
of the fake galaxies m atches the observed i'-drop lu- 
minosity function (e.g., iBouwens et al.l [20061 . For the 
spectral shape, we use synthetic templates (Chariot & 
Bruzual 2010, in preparation) with parameters fixed to 
those which pro vide reasonable fit s for similarly bright 
i'-dropouts (e.g. IStark et "al"1l2009D . Changing these pa- 
rameters to other reasonable values does not affect our 
results. We attach Lya luminosities to each of the galax- 
ies according to an assumed Lya EW distribution (which 
we describe below) and we also account f or Lya fo r est ab - 
sorption using the relations presented in iMeiksinl (|2006f ). 
Finally we derive i' 775 and z§5o-band magnitudes from 
the model SEDs and construct an artificial sample of 
i'-drops which satisfy the color criteria and z-band mag- 
nitude limit. 

We find that the output EW distribution matches the 
input EW distribution of galaxies at the mean redshift 
of the z'-drop population. For example, if we adopt an 
input EW distribution with the form p(WL ya ,o= exp[- 
Wzyafl/Wo] and set Wo=20.0 A, we find that the out- 
put EW distribution is nearly identical to the input dis- 
tribution (Wo=20.1 A). It should be noted that Poisson 
noise (which tends to scatter faint sources toward slightly 
brighter magnitudes) will alter the EW distribution if 
the intrinsic EW distribution is luminosity-dependent, 
as suggested by Paper I. But this effect should occur at 
each redshift and hence should not affect the measured 
redshift evolution in the EW distribution. 

4. RESULTS 

We now derive the EW distribution and Lya fraction 
{XLya) for z ~ 6 galaxies and compare with the lower 
redshift samples of Paper I. We group our i'-drop sam- 
ple into two bins of rest-UV absolute magnitude, taking 
care to apply minor corrections to the observed broad- 
band magnitudes to compensate for the effects of Lya 
emission and IGM Lya forest absorption. For galaxies 
without Lya emission, we correct for Lya forest absorp- 
tion statistically using the redshift distribution predicted 
from the Monte Carlo simulations in §3. 

In Figure 2, we present our observed Lya EW distri- 
bution, with emission lines grouped in 30 A bins. For 
the luminous sub-sample, the distribution, p(M / L yQ .o), 
rises toward lower EW widths, reaching p=12 ± 6.8% in 
the lowest EW bin considered. When compared to the 
EW distribution of luminous sources at 4 < z < 5 (from 
Paper I) , we find that while Lya is marginally more com- 
mon in each EW bin at z ~ 6, the uncertainties are too 
large to distinguish the two distributions. In contrast, 
in the lower luminosity bin, the EW distribution shows 
stronger positive evolution from 4 < z < 5, with Lya 
considerably more prevalent among z ~ 6 LBGs. 

We next compute the LBG Lya fraction by integrat- 
ing the EW distribution above 25 A and 55 A to yield 
the fractions X|^ a and X| 5 aQ . We group galaxies in the 
same two luminosity bins as in the analysis above. In the 
faint subset, we find 54 ± 11% have WL ya ,o > 25 A and 
27 ± 8.0% have W Ly a,o > 55 A. Luminous galaxies ex- 
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Fig. 3. — Predicted rest-frame Lya EW distribution (in bins 
of AW Lya,o = 30 A) for z ~ 7 LBGs based on an extrapola- 
tion of trends from Fig. 2 assuming that the Lya fraction in- 
creases linearly with redshift. Uncertainties are based on statisti- 
cal error in our lower redshift samples. The dashed line indicates 
the limits that could be reached ~ 4 hours of integration with 
Keck/NIRSPEC. Significant deviations below this prediction may 
arise if the IGM is partially neutral. 



hibit Lya emission less frequently, with 20 ± 8.1% and 
7.4±5.0% observed with Lya emission in excess of 25 and 
55 A. Combining our results with those from Paper I, we 
find that the fraction of Lya emitters among the LBG 
population increases with redshift for lower luminosity 
galaxies. Assuming a linear relationship between X^ ya 
and redshift, we find dXf ya /dz = 0.11 ± 0.04. In con- 
trast, less redshift evolution is seen in the larger EW bin 
(dXf ya /dz = 0.018 ±0.036), consistent with the findings 
from Paper I. Similar (albeit noisier) trends are seen in 
the more luminous sub-sample, with the lowest EW bin 
showing the strongest indications of positive evolution 
with redshift. 

This improved determination of the EW distribution 
and Lya fraction for LBGs at z ~ 6 is a necessary step 
toward providing the essential baseline for predicting the 
outcome of spectroscopic campaigns beyond z ~ 7 and 
interpreting any downturn in the Lya fraction that may 
be associated with reionization (see §5). 

5. THE EXPECTED VISIBILITY OF LYa EMISSION IN 
Z > 7 LBGS 

Our new results, taken together with those in Paper I, 
now suggest that ~ 54% of moderately faint (—20.25 < 
Mjjv < —18.75) z ~ 6 LBGs exhibit very strong Lya 
emission. In Paper I, we argued that both the redshift 
and luminosity dependence of the Lya fraction was likely 
due in large part to variations in dust obscuration as 
evidenced by the correlation between Lya EW and the 
rest-frame UV slope, f3. Recent analyses of the colors of 
the z ~ 7 LBGs indicat e that these systems ar e yet bluer 
than those at z ~ 6 (jBouwens et al.l l2010al) . implying 
even less or no dust obscuration. Hence it seems likely 
that the redshift trend in the Lya fraction presented in 
Figure [2] should continue to z ~ 7 and that Lya should 
be readily detectable in sufficiently deep spectroscopic 
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campaigns. 

Given the short cosmic time spanning 6 < z < 7 (~ 170 
Myr), it seems plausible to use the EW distribution and 
Lya fractions presented in Figure [2] to predict the ex- 
pected Lya visibility for sources at z ~ 7, assuming Lya 
flux is not significantly attenuated by neutral hydrogen 
in the IGM. Motivated by the blue z ~ 7 UV slopes dis- 
cussed above, we extrapolate the evolution in XL ya to 
z ~ 7 (Figure 2). For low luminosity sources, this results 
in a small increase in the Lya fraction {&Xf ya = 0.14) 
which we divide into the three EW bins using weights 
set by p(WL ya ,o)- We follow the same procedure for the 
luminous sources. The results, presented in Figure 3, 
suggest a survey of ~ 20-30 galaxies drawn from the now- 
available WFC3/IR target list (e.g.-lMcLure et alj|2010t 
iBouwens et al.ll2010bl:lBunker et alj2010ft would~vield in- 
teresting results. While uncertainty in the observed Lya 
trends and their extrapolation to z ~ 7 obviously af- 
fects our prediction, it seems clear that Lya should be 
common in z ~ 7 samples if the IGM is highly ionized. 
The failure to detect emission in such a sample might 
therefore be a strong indicator of a risin g neutral fraction 
beyon d z ~ 6 as claimed originally bv lKashikawa et all 
(2006) from the luminosity function of LAEs. 

How practical is a search for line emission to the EW 
limit of 20 A discussed above? In terms of an in- 
tegrated line flux FL ya , the EW limit corresponds to 
F^ya ~ 3(7) xl0~ 18 erg cm~ 2 s _1 for z ~ 7 galaxies 
with M[/y=-20(-21) (corresponding to galaxies with ap- 
parent AB magnitudes of J ~ 27 and 26, respectively). 
Such line flux limits are feasible with spectrographs on 
8-10 meter telescopes such as the Near InfraRed SPEC - 
trograph (NIRSPEC) on Keck II ([McLean et al.lfl99l . 
Earlier work with NIRSPEC has reached such a limit 
at 5-cr significance between the atm ospheric sky lines 
in the Y and J-ba nd in ~ 4 hours (|Stark et all 120071 : 
iRichard et a"ill2008[ ). As more multi-object infrared spec- 
trographs become available, it will be feasible to observe 
many z > 7 sources simultaneously, allowing ultra-deep 
exposures of galaxies at least as faint as Muy ~ -19. 

6. CONCLUSIONS 

We present new ultra-deep spectroscopic observations 
of 30 i'-band dropouts in GOODS-N using DEIMOS on 
the Keck II telescope. By adding these spectra to the 
large database of DEIMOS and FORS2 spectra of B, V, 
and i'-band dropouts discussed in Paper I, we demon- 
strate more robustly the rise with redshift over 4 < z < 6 
in the fraction of low luminosity LBGs that show promi- 
nent Lya emission. We also derive a much-improved EW 
distribution of Lya at z ~ 6, the highest redshift where 
the intergalactic medium is known to be highly ionized. 
Motivated by the continued blueward evolution of the 
continuum UV slopes to z ~ 7, we extrapolate the Lya 
fraction trends presented in this paper to redshift 7, and 
predict the likely success rate of recovering Lya emis- 
sion in the new population of z > 7 sources located with 
HST. Our results suggest line emission should be readily 
detected given adequate observing time and that quan- 
titative results would therefore test for the presence of 
neutral gas associated with the end of cosmic reioniza- 
tion near z ~ 7. 
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